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Abstract 

The existence of two sterile neutrinos Ug and z/' (blind to all Standard-Model inter¬ 
actions) is shown to be implied by a model of fermion ’’texture” that we develop since 
some time. They may mix nearly maximally with two of three conventional neutrinos, 
say Ue and thus leading to neutrino oscillations, say z/g —>■ Vg and z/^ —>• z/(, with nearly 
maximal amplitudes. Then, they can be responsible for the observed deficits of solar izfs 
and atmospheric z/^’s, respectively, but by themselves do not help to explain the LSND 
results for oscillations. On the other hand, they are consistent with the CHOOZ 

negative result. At the moment, the experiment cannot decide, whether the deficit of 
atmospheric z/^’s, confirmed by the recent Super-Kamiokande findings, has to be related 
to the oscillations v^^ —>■ or z/^ —>■ z/(. In the last Section of the paper, a new notion 
of ’’non-Abelian spin-1/2 fermions” is presented in the context of a composite option for 
fermion families. 

PACS numbers: 12.15.Ff , 12.90.-|-b , 14.60.Gh 


August 1998 


Dedicated to Jan Lopuszahski in honour of his 75th birthday. 



1.Introduction 


The hypothetic sterile neutrinos, by dehnition interacting only gravitationally, are 
blind to all Standard-Model interactions, in contrast to the conventional neutrinos (or, 
rather, their lefthanded parts) which participate hrst of all in the weak sector of Standard- 
Model interactions. Such Standard-Model-inactive fermions are invoked from time to 
time by theorists, who want to explain [1] through neutrino oscillations not only the 
observed dehcits of solar and atmospheric neutrinos, but also the results of LSND exper¬ 
iment. The sterile neutrinos may also form a Standard-Model-inactive fraction of the 
dark matter. 

In the present paper, we demonstrate how two different sterile neutrinos are implied 
by a model of fermion ’’texture” [2,3] that we develop since some time. As shown previ¬ 
ously, this model justihes [2] the existence of three and only three families of conventional 
leptons and quarks (z/e , e~ , u , d), (z/^ , , c, s), {ur , t~ , t, b) and, moreover, describes 

reasonably [3] the masses and mixing parameters of quarks and charged leptons, making 
also some useful suggestions as to neutrinos. Note that in this model all neutrinos are 
Dirac particles having both lefthanded and righthanded parts. 

In order to make our presentation fairly comprehensible, we will hrst recapitulate 
briehy the basic features of the model in its part concerning the existence of fundamental- 
particle families [2]. Then, we shall discuss the existence of two sterile neutrinos and the 
related neutrino oscillations. 


2. Dirac’s generalized square root 

The starting point of our model is the conjecture that a//kinds of matter’s fundamental 
particles existing in Nature can be deduced from Dirac’s square-root procedure 
T ■ p. 

As is easy to observe, this procedure leads in general to the sequence N = 1, 2, 3,... 
of different (generally reducible) representations 


of the Dirac algebra 


T^ = 


1 N 

Tv S'’-" 


( 1 ) 


1 





( 2 ) 


constructed with the use of the sequence iV = 1, 2, 3,... of Clifford algebras 

= (i,j = l,2,...,]V). (3) 

Then, the sequence iV = 1, 2, 3,... of Dirac-type equations follows, 

{T-\p-gA{x)]-M}i:{x) = 0, (4) 

where gT ■ A{x) may symbolize the minimal coupling of ip{x) to the Standard-Model 
gauge helds A^{x) including all SU{3) x SUl{2) x U{1) coupling matrices: A’s, r’s, Y 
and r® = 

In Eqs. (4) the matrices (1) can be presented in the reduced forms 

= 1 0 • • • 0 1 (5) 

(iV—l)times 

with 7 ^ and 1 denoting the usual 4x4 Dirac matrices. Then, the Dirac-type equations 
(4) can be rewritten as 


{7 ■ [p - gA{x)] - ^/ 3 ia 2 -aiv = 0 ( 6 ) 

with 'ip{x) = {'ipaia 2 -aN{x)), where ai, a 2 ,---, ctat stand for N Dirac bispinor indices: 
ttj = 1, 2, 3,4 for i = 1,2,..., N. Here, the chiral representations are used to dehne all 
ttj {i = 1,2,..., N). This means that ai = 1,2, 3,4 correspond to four different pairs 
(1,1), (1,-1), (-1,1), (-1,-1) of eigenvalues of the matrices 


Tf = iT^T^r^^r)* , = ir)T)T 


''0'pl'p2'p3 


3 _ 


''5'p0'p3 
i ’ 


( 7 ) 


simultaneously diagonal for all i, which choice is allowed because all T^ and Ef commute 
both for equal and different i. The Tf matrices (i = 1, 2,..., N) appearing in Eqs. (7) are 
dehned as N (properly normalized) Jacobi combinations of yf matrices {i = 1,2,..., N), 
where in particular Tj* = T^ is given as in Eq. (1). Then, {Tf, Tj} = 26ijg^'^ {i,j = 


1,2,..., N) due to Eqs. (3), and also {Tf , TM = 0, but 


■p5 "pS 
i ’ i 


= 0, where particularly 


2 



= r®. Note that in the one-body Dirac-type equations (4) there appear only the 
”centre-of-mass” matrices, while all ’’relative” matrices r 2 ,...,r^ are absent. In 
spite of this, all ai,a 2 , ■ ■ ■ ,q:n are present in Eqs. (6): both the ” centre-of-mass” Dirac 
bispinor index ai as well as the ’’relative” Dirac bispinor indices 02 , • • • ,C(n, the latter 
are decoupled, however, even in the presence of Standard-Model coupling gTi ■ A[x). 

For iV = 1 Eq. (6) is obviously the usual Dirac equation, for iV = 2 it is known as 
the Dirac form [4] of the Kahler equation [5], whilst for N = 3,4, 5,... we obtain new 
Dirac-type equations [2], 

If the Standard-Model coupling gV ■ A{x) is really present in Eqs. (6), then the Dirac 
bispinor index ai, which is the only ai affected by the gauge helds A^{x), is distinguished 
by its correlation with the set of all diagonal Standard-Model charges ascribed to any 
particle of the helds t/’ai a 2 ...ajv(^) label / of this set is here suppressed). The remaining 
Dirac bispinor indices 0 : 2 ,..., aw are all decoupled and so, physically unobservable in the 
gauge helds A^(a;). It is natural to conjecture that they are physically undistinguishable 
and, therefore, are formal objects obeying Fermi statistics along with Pauli principle. This 
implies that t/’oi a 2 ...a]v(^) antisymmetric with respect to a 2 ,..., aw- 

The above conjecture, together with the probabilistic interpretation of wave functions 
t/’oi 02 ...aiv(^) requirement of their relativistic covariance applied to all bispinor 

indices ai, a 2 ,..., aw, leads to the conclusion that there are three (and only three) families 
N = 1, 3, 5 of leptons and quarks [2], and two (and only two) families N = 2, 4 of some, 
not yet observed, fundamental scalars [6]. They correspond to the wave functions 


and 




T {c y \ t(’oi0203 

4 ^ ^ 0i.20i.Z 

£q ,2 Qg 04 Qg'^ai 02 as a4 0^5 


t/’oi 12 — ll’ai 34 , 


t/’ai 1234 


( 8 ) 



2 V 2 

1 


(C ^aia2 = (^12 - ^21) = 4= (^34 “ ^43) , 

\ /aia2 4/2 V2 


4^0 


'a.\a. 2 <y.z^A: r q;iQ; 2<^3<^4 




(t/’l234 ~ t/’2134 + t/’3412 ~ t/’4312) , 


(9) 
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respectively. Each of these wave functions carries the (here suppressed) Standard-Model 
label / = u , e, u, d denoting four sorts of fundamental particles corresponding to the 
signature of conventional neutrinos u and charged leptons e as well as up quarks u and 
down quarks d, all four following from the Standard Model (though the existence of three 
and two fundamental-particle families does not follow from it). In the case of fundamental 
fermions, the three families are, of course, (z/g , e“ , u , d), , fi~ , c, s) , {ur , t~ , t, b) , 

while in the case of fundamental scalars one of (a priori) possible options may be that 
the two families correspond to the first and second fermion family [ 6 ]. 

Now, in contrast, if the Standard-Model coupling gT ■ A{x) is absent from Eqs. 
(6), then only physically undistinguishable he., antisymmetric bispinor indices ai {i = 
1,2,..., N) can appear at the wave functions (^)- ^his case, the argument 

similar to the used before shows that on the fundamental level there are two (and only 
two) Standard-Model-inactive spin-1/2 fermions iV = 1,3 corresponding to the wave 
functions 





Oil 0.2 


-020^04,0^ YO^O^O^ 




( 10 ) 


(with no suppressed / label). They can be identihed with two sterile neutrinos denoted in 
this paper by Ug and z/(, respectively. Analogically, on the fundamental level there should 
exist also two (and only two) Standard-Model-inactive spin-0 bosons = 2 ,4 that may 
be called sterile scalars, 0 *^^^ and 0 *^^^ (with no suppressed / label). 


3. Neutrino oscillations involving Ug and Ug 


Let us conjecture tentatively that the sterile neutrinos Ug and z/( are compelled to 
mix nearly maximally with the conventional neutrinos z/g and z/^, respectively, in order 
to form four related neutrino mass states z/i or z /4 and z /2 or z/ 5 . Other neutrino mixings 
are assumed not to appear at all or to be negligible. In particular, the third conventional 
neutrino is left not mixed and so, z /3 = z/,- is a neutrino mass state. Evidently, the 
mixings of z/g with Ug and z/^ with z/( would be forbidden, if the electroweak SUl(2) x 17(1) 
symmetry were not spontaneously broken. Thus, we can say that neutrino oscillations, 
being a consequence of these mixings, are caused in fact by the spontaneous breaking of 


4 



electroweak symmetry (if, of course, sterile neutrinos exist). 

Under the above conjecture, the unitary transformation vj = Y,a Via^a between neu¬ 
trino mass states z// = z/i, z/ 2 , r's, 1 ^ 4 , 1^5 and neutrino flavor states z/q, = z/g, z/^, i^r, ^s, K 
is given as 


Z/i = UiiZ/e + U14Z/S 
z/2 = U22Z/^ + V'25Z/( 
^3 


Z/4 = U4iZ/e H- U44Z/S , 
^5 = + U55Z/( , 

Z/^ , 


where the nonzero coefficients are 


Uii = U 44 = 
U 22 = U 55 = 


1 

1 

v^IT^ 

U 33 


Ul4 = 

1 

* 

Y 

^l + Y^ 

U25 = 

—V* — - 

*^52 — 

- ^ e-' 

V1+X2 


( 11 ) 


( 12 ) 


(in Eqs. (11) and (12), for notation convenience, we write Vjj in place of Vja, where 
/, J = 1, 2, 3,4, 5). The magnitudes of these coefficients are determined by the parameters 


Y 


Mn - _ M 44 - 

IM14I |Tfi 4 | 


X 


M 22 - 

IM 25 I 


M55 - 771, 

IM 25 I 


(13) 


involving neutrino masses 


m 


iff44 




T 


— M 4.4 


+ \M, 


141 


m 


M22 + M55 


142,145 


M 22 — M; 


55 


+ \M, 


251 


(14) 


On the other hand rriu^ = M 33 . Here, {Mu) {I,J = 1,2, 3,4, 5) is a 5 x 5 neutrino 
mass matrix with M 14 = = |Mi 4 | exp 799 and M 25 = = I "^ 25 1 exp 799 ' as the only 

off-diagonal elements. Then, (Vj-j) {I,J = 1,2, 3,4, 5) is a 5 x 5 lepton counterpart of 
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the familiar Cabibbo—Kobayashi—Maskawa matrix for quarks, where now Vu = 
and V 25 = —V ^2 ^he only nonzero off-diagonal elements. 

Some (here neglected) small corrections to the neutrino mixings (11) may be caused by 
possible small deviations of the charged-lepton mass matrix from a diagonal form [3]. In 
fact, these deviations produce small deviations of the related diagonalizing unitary matrix 
from the unit matrix. In turn, such a charged-lepton diagonalizing matrix contributes 
multiplicatively to the lepton Cabibbo—Kobayashi—Maskawa matrix [3], changing a little 
its leading form ( 12 ) (in particular, nearly all zero elements of its leading form become 
nonzero but small). 

Now, making use of Eqs. (12), we can calculate the probabilities of neutrino oscil¬ 
lations z/e —^ i^s and z/( (in the vacuum) from the general formula 

/ ArT?2 \ 

P{i^a = IWka(t))P = ^ I4^K;„VJ^Kxaexp 

where (on the rhs of Eq. (15), for notation convenience, we will 

replace a, P hj I, J = 1,2, 3,4, 5). Here, z/a(0) = z/q,, {z/^| = (0|z/^, (zz/jIz/q) = and, as 
usual, t/\p\ = L/E (c = 1 = h), what is equal to 4 x 1.2663L/E if , L and E are 

measured in eV^, m and MeV, respectively {L is, of course, the source-detector distance). 
Further on, we will denote 

xlk = 1.2663^^^^!^ (16) 

tj 

and use the identity cos 2 xlk = 1 — 2 sin^ xlk- 

In such a way, we derive the following formulae for probabilities of neutrino oscillations 
z/g —z/s and z/^ —z/( (in the vacuum): 


^2 

P ^ i^s) = sin^ a;4i , 

= ^ (1 + X 2)2 , ( 17 ) 

while all other P (z/^ —> z/^) with a 7 ^ /? vanish [except, of course, for P (z/^ —z/g) = 
P (z/g — Ug) and P (z/( —z/^) = P (z/^ — u'P)]. Thus, the neutrino-oscillation formulae (in 
the vacuum) for survival probabilities of z/g and z/^ are 
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Y 2 

J^2 

P ^l^p,) = 1 - 4 ^^ ^^ 2)2 ^52 • (18) 

In the case of solar neutrinos, the observed dehcit of u^s can be explained through 
the neutrino oscillations (in the vacuum), when using the two-flavor formula for survival 
probability of z/g. 


PU 


with the parameters [7] 


1 sin 20 sQ[a^j. sin 





(19) 


Sin2 20soiar ~ 0.65 to 1 , ~ (5 to 8 ) x 10"^^ eV^ . (20) 

These give the so called vacuum fit, in contrast to two other known fits based on the 
resonant MSW mechanism [ 8 ] in the Sun matter. In our model of neutrino oscillations 
(where Vg oscillations are responsible for the dehcit of solar z/g’s), this ht leads to 


4y2 

--—^ ~ 0.65 to 1 , Am?, ~ (5 to 8 ) x 10"^^ eV^ (21) 

(1 + r ^)2 

(as m^^ > m^J. Hence, Y ~ 0.507 to 1 and so, we get a large mixing of z/g with z/^: 
Vii = V 44 ~ 0.892 to 1/v^ and V 74 = —~ —(0.452 to l/\/2)expi(f (the phase (f 
remains not determined). 

In the case of atmospheric neutrinos, the recent hndings of the Super-Kamiokande 
experiment [9] show that the observed dehcit of z/^’s can be explained also through the 
neutrino oscillations (in the vacuum), when making use of the two-havor formula for 
survival probability of z/^. 


P (z2^ ^ Z2^) 


= 1 — sin^ 26'atm sin^ 


B 


with the parameters 


( 22 ) 


sin2 20atm ~ 0.82 to 1 , ~ (0.5 to 6 ) x 10-^ eV^ . (23) 
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In our model of neutrino oscillations (where z/( oscillations are responsible for the 

dehcit of atmospheric z/^’s), this implies 

4X^ 

~ to 1 , Am 52 ~ (0.5 to 6) x 10“^ eV^ (24) 

(as > ml^). Hence, X ~ 0.636 to 1 and thus, we obtain a large mixing of z/^ with z/(: 
^22 = ^55 ~ 0.844 to 1/v^ and V 25 = —h ^2 ~ —(0.537 to l/v^)expi99' (the phase (f' 
remains not determined). 

On the other hand, the CHOOZ experiment [10] found no evidence for neutrino- 
oscillation modes of z?e in a parameter region overlapping the range (23) of sin^ 26'atm and 
what shows that within this parameter range there are no mixings of i/g with z/^, 
z/t-, i/g, z/(. In particular for z/^, this is consistent with the assumed dominance of mixing 
between z/^ and z/( over mixing between z/^ and z/g within the range (23) of sin^ 20atm and 
Am^^m the moment, however, it cannot be decided experimentally [9], whether the 
mixing of z/^ with z/( or the here neglected mixing of z/^ with i/j- is responsible for the 
dehcit of atmospheric z/^’s). For i/g, this requires that the assumed strong mixing of z/g 
with i/g must correspond to parameters sin^ 26'soiar and Arrig^i^j. belonging to a range very 
different from (23) [in fact, they can he in the range (20)]. Finally for z/,-, the lack of 
mixing between and z/g is one of necessary and sufficient conditions for the assumed 
identity 1/3 = Vr (another is the lack of mixing between and z/^, if this really is true). 

Of course, the sterile neutrinos Vg and z/( by themselves cannot help to explain the 
results of LSND experiment [11] which gave evidence for z/^ —z/g and z/^ —z/g oscillations 
corresponding to a much larger Am^gNo than both and These oscillations, 

if eventually conhrmed, would require a considerable mixing of z/^ with z/g, corresponding 
to parameters sin^26*LSND and Am^gNo a range very different from (23). This 

mixing should be stronger than that induced by the (mentioned before) nondiagonal 
charged-lepton corrections appearing in our model of fermion ’’texture” [3]. 

The comparison of mass squared differences Amh and Am |2 as estimated in Eqs. 
(21) and (24) suggests that and are possibly much smaller than and 
(alternatively, and ml^ may be much more degenerate than ml^ and ml^). 


4. Outlook: Non—Abelian spin—1/2 fermions 



When the Dirac-type equations (4) are considered, one may ask a (perhaps) profound 
question, as to whether these one-body equations could be understood physically as 
point-like limits of some iV-body equations for tight bound states of N spin-1/2 preons 
with equal masses. If it was so, the four-positions of such subelementary constituents 
(called here preons, as usual) should tend practically (within the bound states) to their 
centre-of-mass four-position, 

N N 

Xi = X + 5xi ^ X , X = ^ ^= 0 , (25) 

i=i i=i 

while then Spi, defined by their four-momenta 

N N 

Pi = P + 5pi , P = J2Pi , = ^ ^ (26) 

i=l i=l 

should vanish in action on the wave functions [here, Xi = {ti , Xi), 5xi = ((5fj, 5xi) and 
X=(f,X)]. 

Of course, the physical mechanism for realization of such practically point-like limits in 
A^-body systems would be provided by an unknown, very strong and shortrange attraction 
between their N constituents (described, for convenience, in the equal-time formalism, 
where 5ti = 0 and 5p^ vanish in action on the wave functions). The (necessarily) non- 
Standard-Model nature of such an attraction would be certainly the most obscure aspect 
of the compound option for the Dirac-type equations (4). 

Let us denote by Pi and Xi {i = 1,2, ...,iV), with Pi = P and Xi = X, the 
(properly normalized) Jacobi combinations of four-momenta pi and four-positions Xi 
{i = 1, 2,..., iV), respectively, for N particles. Then, 

pr.-q] (27) 

Making use of this notation, we can write the identities 

( 7 i ■ Pi - nii) = —= Y (r* ■ Pi - \fNmi) , (28) 

i=l vA/ ^ 

where Tf (i = 1,2,...,A^), with T^ = T^, stand for the (properly normalized) Jacobi 
combinations of yf matrices (i = 1, 2,..., N) for N particles [the Tf matrices were already 
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introduced in Eqs. (7), though only in reference to the one-body Dirac-type equations(4)]. 
Then, 


{rf, r;} = 25,,{t,j = 1,2,..., iv), (29) 

as follows from Eqs. (3). Here, in particular, T^ = T^ is given as in Eq. (1). 

In this notation, the natural candidates for the hypothetic iV-body equations would 
be 


|ri ■ [Pi - gA{X,)] + f; r, ■ P, - I ^(Xi, X 2 ,..., Xiv) = 0 , (30) 

where /(X 2 ,... ,Xiv) would symbolize the unknown non-Standard-Model attraction be¬ 
tween N constituents. In Eqs. (30), the Standard-Model gauge fields A^[x) are coupled 
to the hypothetic iV-body systems at four-points describing their centre-of-mass four- 
positions Xi = X. This is approximately true, when A^{X + Sxi) are only weakly 
dependent on 6 xi. 

In the point-like limits, where the relative four-positions X 2 ,..., Xjsf [i.e., also all 
5xi) tend to zero and then the relative four-momenta P 2 ,..., Pn {i.e., also all 6 pi) vanish 
in action on the wave functions, Eqs. (30) are really reduced to the Dirac-type equations 
(4) withp = P = Pi, X = X = Xi, T = Ti and M = \/N {Nm + Ixi^o) {rrii = m). Note 
that M grows with N faster than linearly. 

In the equal-time formalism, where the relative times t2, ■ ■ ■ An (he., also all Sti) are 
zero and the relative energies P 2 , ■ ■ ■, P^ (he., also all Sp^) vanish in action on the wave 
functions, Eqs. (30) assume the forms 


N 


PiA{Xi,X2 , . . . ,XN,t) — \TlTi-\Pi-gA{Xi,t) \ +gA\Xi,t) + J2T%-P, 


i=2 


+T«(VN'£ m, + j I V>(Vi, W, ..., i) . 


( 31 ) 


where P° = P° = id/dt and Ixo=o = h(X 2 , • • • ,Xx)- 

In the point-like limits, Eqs. (31) are reduced to the equations 
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t) = |r°f ■ p — gA{x,t) + gA^{x,t) + (32) 

with p = P = Pi, X = X = Xi, r = ri and M = \/N{Nm + Ixi^o) = m). Of 
course, p^ = P^ = id/dt and Ixi^o stands for a reasonably defined point-like limit of J. 

Note that the eigenvalues (P°)kin of the kinetic part of the hamiltonian appearing on 
the rhs of the state equation (31) get for any N the form ± P^ + ... + P^ + iV(iVm)^ 

= ±y/N [pf + • • • + Pat + iV(iVm)^]^^^, as if our iV-body system were a single Dirac parti¬ 
cle with the mass Nm in a (3iV-|-l)-dimensional spacetime (notice, however, the additional 
factor V^)- 

A fundamental feature of Eqs. (30) is that, via Ff {i = 1,2,..., N), they contain N 
Dirac nonconventional yf matrices (f = 1, 2,..., N) which do not commute for different 
particles, in contrast to Dirac conventional gammas commuting for different particles [in 
fact, the nonconventional yf (i = 1, 2,..., N) anticommute for different i, as is seen from 
Eqs. (3)]. The spin-1/2 fermions i = 1,2,...,N described within an iV-body system 
with the use of such nonconventional yf matrices {i = 1,2,..., N), anticommuting for 
different particles, might be called non-Abelian spin-l/2 fermions [12], In contrast, in the 
familiar case of Dirac conventional gammas, commuting for different particles, one could 
use the term Abelian spin- 1/2 fermions. 

Now, let us observe that the form of spin tensors for spin-1/2 fermions f = 1, 2,..., iV 
is identical in the non-Abelian and Abelian case: 


fiu _ ^ r /i. i/i _ 

~ 2 1 for fi = k, u = I ’ 

where = y/y[ and a™ = y/y/y™ = yfc^r with y/ = i^^'ll if if ■ In fact, for each i the 
components satisfy in both cases the usual Lorentz-group commutation relations, 
while for different i they commute in both cases as being bilinear in yf . Also if commute 
for different i in both cases. So, the total generators of Lorentz group for a system of N 
spin-1/2 fermions are given in both cases by the operators 


for /i = 0 , z/ = / 




Let us note, by the way, the following identity valid for the total spin tensor in both 
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( 35 ) 


N N 

i:<rr = Esr. 

i=l i=l 

where 


Ef^-irr.rci^ 


iA\ for /i = 0, z/ = / 

eklmj^m fQj, jji = k, V = I 


(36) 


with A\ = r°r( and Ef = rfr°r™ = TfA^. Evidently, Sf = E^^ with E'^^ = | [E^, T^] is 
the centre-of-mass spin tensor for the system of N spin-1/2 fermions, while E^^,..., E)^^ 
are its relative spin tensors. All spin tensors E/*^, being bilinear in F/, commute for 
different i in both cases. 

Thus, in this Section, we can draw the important conclusion that for a system of 
N spin-1/2 fermions the Lorentz-group commutation relations get two (and only two) 
realizations: either with the use of Dirac conventional gammas commuting for different 
particles, or with the use of Dirac nonconventional gammas anticommuting for different 
particles. Such an intriguing statement seems to support the logical consistency and unex¬ 
pected naturalness of the notion of non-Abelian spin-1/2 fermions. They may provide an 
unconventional alternative for familiar Abelian spin-1/2 fermions in the potential struc¬ 
ture of particle theory. In this Section, their role as hypothetic preons was underlined. 

Finally, we should like to emphasize some unconventional features of the quantization 
procedure which would work in the case of non-Abelian spin-1/2 fermions. It is not dif- 
hcult to observe that in the case of spin-1/2 fermions, only the Fock-space states related 
to Dirac conventional gammas (commuting for different particles) can be constructed by 
means of the familiar second-quantization procedure based on Fermi creation and anni¬ 
hilation operators for single particles. This is so, because the repetition of single-particle 
creation operators can lead to Fock-space states of particles with commuting Dirac gam¬ 
mas only. In order to construct the Fock-space states related to Dirac nonconventional 
gammas (anticommuting for different particles), new Fermi or Bose operators creating and 
annihilating at once whole A^-particle configurations with odd or even N = 1,2,3,..., 
respectively, must be introduced. Of course, these N particles are then non-Abelian 
spin-1/2 fermions. Such a new procedure might be called the third quantization [12]. 
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